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Abstract 

This paper addresses some fundamental issues in nonconvex analysis. By using pure 
complementary energy principle proposed by the author, a class of fully nonlinear par¬ 
tial differential equations in nonlinear elasticity is able to converted a unified algebraic 
equation, a complete set of analytical solutions are obtained for 3-D finite deformation 
problems governed by generalized neo-Hookean model. Both global and local extremal 
solutions to the nonconvex variational problem are identified by a triality theory. Con¬ 
nection between challenges in nonlinear analysis and NP-harcl problems in computational 
science is revealed. Results show that Legendre-Hadamard condition can only guarantee 
ellipticity for generalized convex problems. For nonconvex systems, the ellipticity depends 
not only on the stored energy, but also on the external force field. Uniqueness is proved 
based on a generalized quasiconvexity and a generalized ellipticity condition. Application 
is illustrated for nonconvex logarithm stored energy. 
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1 Nonconvex Variational Problem and Challenges 

Minimum total potential energy principle in nonlinear elasticity has always presented funda¬ 
mental challenging problems not only in continuum mechanics, but also in nonlinear analysis 
and computational sciences. This paper intends to solve, under certain conditions, the fol¬ 
lowing minimum potential variational problem {(V) for short): 

(V) : min jn( X ) = jf W(V X )dB - jf * • tdS| A G Y c j , (1) 

where the unknown deformation y(x) = {Xi( x j)} 6 Y a is a vector-valued mapping B C M 3 —> 
uj C R 3 from a given material particle x = {xi} £ B in the undeformed body to a position 
vector in the deformed configuration uj. The body is fixed on the boundary S x C dB, while 
on the remaining boundary St = S x n dB, the body is subjected to a given surface traction 
t(x). In this paper, we let X a as a geometrically admissible space defined by 

Ya = {A G W 1 ’ 1 (B; M 3 ) | *(x) = 0 Vx g S x } (2) 

where W 1,1 is the standard notation for Sobolev space, i.e. a function space in which both a 
and its weak derivative Va have a finite L 1 (S) norm. For homogeneous hyperelastic body, 
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the strain energy W(F) is assumed to be C l on its domain T c C M 3x3 , in which certain 
necessary constitutive constraints are included, such as 


detF > 0, W(F) > 0 VF £ T c , W( F) -> oo as ||F|| —> oo. 
Thus, the kinetically admissible space in (V) is simply defined by 

x c = {x £ X a \ V* £ -F c } 


(3) 

(4) 


which is essentially nonconvex due to nonlinear constraints such as det(V%) > 0. Also, the 
stored energy W (F) is in general nonconvex in order to model real-world problems such as 
post-buckling and phase transitions, etc. Therefore, the nonconvex variational problem (V) 
has usually multiple local optimal solutions. 

Let X b C X c be a subspace with two additional conditions 


X b = {xe X c \ X G C 2 (B-R 3 ), W(Fix)) G C 2 (P c ; M)}. (5) 

If dB is sufficiently regular, the criticality condition 5II(x; 5x) = 0 £ X b leads to a 

nonlinear boundary-value problem 


(BVP) : 


f -V • cr(Vx) = 0 in B, 

\ N • it(Vx) = t on S t , X = 0 on S x 


( 6 ) 


where, N £ M 3 is a unit vector normal to dB, and <r(F) is the first Piola-Kirchhoff stress 
(force per unit undeformed area), defined by 


<r = VW(F), 


or dij = 


dW(F) 

~dF~ 


i,j = 1,2,3. 


(7) 


Remark 1 (Nonconvexity, Multi-Solutions, and NP-Hard Problems) 

The stored energy W(F) in nonlinear elasticity is generally nonconvex. It turns out that 
the fully nonlinear (BVP) could have multiple solutions {% fc (x)} £ X c C M°° at each ma¬ 
terial point x £ B s C B. As long as the continuous domain B s / 0, this solution set 
{x fc (x)}(fc = 1,..., K) can form infinitely many ( K °°) solutions to (BVP) even B Cl. It is 
impossible to use traditional convexity and ellipticity conditions to identify global minimizer 
among all these local solutions. Gao and Ogden discovered in [10j that for certain given 
external force field, both global and local extremum solutions are nonsmooth and can’t be 
obtained by Newton-type numerical methods. Therefore, Problem ( V) is much more difficult 
than (BVP). In computational mechanics, any direct numerical method for solving ( V ) will 
lead to a nonconvex minimization problem in M n , which could have K n local solutions. Due 
to the lack of global optimality condition, it is fundamentally difficult to solve nonconvex 
minimization problems by traditional methods within polynomial time. Therefore, in com¬ 
putational sciences most nonconvex minimization problems are considered to be NP-hard 
(Non-deterministic Polynomial-time hard) (Tlj. 

Direct methods for solving nonconvex variational problems in finite elasticity have been 
studies extensively during the last fifty years and many generalized convexities, such as poly-, 
quasi- and rank-one convexities, have been proposed. For a given function W : T c —> M, the 
following statements are well-known (see |18| Q 

convex =>• poly-convex => quasi-convex =>- rank-one convex. 

x It was proved recently that rank-one convexity also implies polyconvexity for isotropic, objective and 
isochoric elastic energies in the two-dimensional case B51- 
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Although the generalized convexities have been well-studied for general function W(F) on 
matrix space W nxn , these mathematical concepts provide only necessary conditions for local 
minimal solutions, and can’t be applied to general (nonconvex) finite deformation problems. 
In reality, the stored energy W{ F) must be nonconvex in order to model real-world phenom¬ 
ena. Strictly speaking, due to certain necessary constitutive constraints such as det F > 0 and 
objectivity condition etc, even the domain T c is not convex, therefore, it is not appropriate to 
discuss convexity of the stored energy W (F) in general nonlinear elasticity. How to identify 
global optimal solution has been a fundamental challenging problem in nonconvex analysis 
and computational science. ■ 

Remark 2 (Canonical Duality, Gap Function, and Global Extremality) 

The objectivity is a necessary constraint for any hyper-elastic model. A real-valued function 
W : T c -> M is objective iff there exists a function V (C) such that !F(F) = V (F r F) VF e T c 
(see CD- By the fact that the right Cauchy-Green tensor C is an objective measure on a 
convex domain £ a = {C e M 3x3 | C = C T , C >- 0}, it is possible and natural to discuss 
the convexity of F(C). A real-valued function V : £ a —> M is called canonical if the duality 
relation = VF(£) : £ a —> £* a is one-to-one and onto [5]. The canonical duality is necessary 
for modeling natural phenomena, which lays a foundation for the canonical duality theory 
|5J. This theory was developed from Gao and Strang’s original work in 1989 [ 12j for general 
nonconvex/nonsmooth variational problems in finite deformation theory. The key idea of this 
theory is assuming the existence of a geometrically admissible (objective) measure £, = A(F) 
and a canonical function F(£) such that the following canonical transformation holds 

£ = A(F) : T a -> £ a =» W(F) = V(A(F)). (8) 

Gao and Strang discovered that the directional derivative A<(F) = <5A(F) is adjoined with 
the equilibrium operator, while its complementary operator A C (F) = A(F) — A*(F)F leads 
to a so-called complementary gap function , which recovers duality gaps in traditional duality 
theories and provides a sufficient condition for identifying both global and local extremal 
solutions mm- ■ 

The canonical duality theory has been applied for solving a large class of nonconvex, 
nonsmooth, discrete problems in multidisciplinary fields of nonlinear analysis, nonconvex 
mechanics, global optimization, and computational sciences, etc. A comprehensive review 
is given recently in m- The main goal of this paper is to show author’s recent analytical 
solutions [7] for general anti-plane shear problems can be easily generalized for solving finite 
deformation problems governed by generalized neo-Hookean materials. Some insightful results 
are obtained on generalized convexity and ellipticity in nonlinear analysis. 

2 Complete Solutions to Generalized Neo-Hookean Material 

By the fact that the right Cauchy-Green strain C = F 7 F is an objective tensor, its three 
principal invariants 

/ 1 (C) = trC, J 2 ( C) = ^[(trC) 2 -tr(C 2 )], J 3 (C) = detC (9) 

are also objective functions of F. Clearly, for isochoric deformations we have / 3 (C) = 1. 
The elastic body is said to be generalized neo-Hookean material if the stored energy depends 
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only on I\, i.e. there exists a function V(I\) such that VF(F) = V(/i(C(F))). Since I\ = 
tr(F 7 F) > 0 VF E P c , the domain of V(I\) is a convex (positive) cone 

Sa = £ LP(B) I £(x) > 0 Vx E B}, (10) 

it is possible to discuss the convexity of V (I\) on £ a . Furthermore, we assume that V (I\ ) 
is a C 2 (S a ) canonical function. Then the canonical transformation (|8j for the generalized 
neo-Hookean model is 

£ = A(F) = tr(F T F) : F c ^£ a , W(F) = F(£(F)). (11) 

For a given external force t(x) on St, we introduce a statically admissible space 

% = {TeW 1 ’ 1 (6;l 3x3 ) | V-T = 0 in B, N • T = t onS t }. (12) 

Thus for any given T E T a , the primal problem (V) for the generalized neo-Hookean material 
can be written in following canonical form 

(F)t : min |n T (V X ) = J G(V X ) dB | V* E T c j , (13) 

where X c = {^E X a \ A(V^) £ £ a } and the integrand G : T a —> M is defined by 

G{F) = F(A(F)) -tr(F r T). (14) 

By the fact that det F > 0 is not a variational constraint and the certain constitutive con¬ 
straints, such as coercivity and objectivity, have been naturally relaxed by the canonical 
transformation, the domain of G(F) is simply T a = M 3x3 . 

Let SO(3) = {R E M 3x3 | R t = R 1 , detR = l}and 

n = {R(x) E L 1 [H,M 3x3 ]| R(x) E SO(3) Vx E B}. (15) 

Theorem 1 For any given T E T a , if X £ is a stationary solution to (V)t, then it is also 
a stationary solution to (' P ). 

For any given rotation field R(x) £ 1Z such that R 7 T E T a , then Ht(F) = Ht(RF). 

For any uniform rotation R E SO(3) such that R 7 T E T a , if X is a stationary solution to 
(V), then R% is also a stationary solution to (V). 


Proof. For any given T E Ta, the stationary condition for the canonical variational problem 
(V)t leads to the following canonical boundary value problem 


(BVP ) T : 


f V • (2CV*) = V • T = 0 in B, 

\ N • (2CV*) = N • T = t on S t , X = 0 on S x 


which is identical to (BVP) since 

<t = V1F(F) = ^P^ = 2CF, C = W(0- 

By the objectivity of ^ = A(F) = A(RF) VR(x) E 1Z and the fact that 

[ tr[(RV X ) T T]dH= [ tr[(Vx) T (R T T)]d£= [ x -tdS VR t T E T a , 
JB Jb JSt 


(16) 
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we have IIx(F) = Llx(RF) VR(x) £ 1Z. Particularly, for any uniform R £ S0(3) such that 
R r T £ 7;, we have II(x) = n T (RF(x)). □ 


Theorem [T] is important for understanding the canonical duality theory. 

By the canonical assumption on V(£), the duality relation £ = W(£) : £ a —> £* a is invert¬ 
ible. The complementary energy can be defined uniquely by the Legendre transformation 

F*(C) = {£C-F(OIC = VF(£)}. (17) 

Clearly, the function V : £ a — > M is canonical if and only if the following canonical duality 
relations hold on £ a x £* 

C = W (0 4 * ? = vu*(C) 44 v(0 + v*(0 = ac (18) 

Using V(£) = — U*(C)> the nonconvex function G(F) can be written as the standard Gao 

and Strang total complementary function E : X a X £* — > M 

3(X,C)= [ [A(V X )C-F*(C)-tr((Vx) T T)]dR (19) 

Jb 

Let S a C5*bea canonical dual feasible space defined by 

S a = {(e£* a \(- 1 T 2 eL\B)}. ( 20 ) 


Then for a given £ £ S a , the canonical dual function can be obtained by the canonical dual 
transformation: 

R d (C) = sta{H(x, C)| X&X a }= t G d (()dB , (21) 

JB 

where the notation sta{S(x, 01 X £ X a } stands for finding (partial) stationary point x £ T 0 
of E(x, 0 for a given £ £ S a , and 

G d (0 = -V*(C) - |CV, T 2 = tr(T T T). (22) 

Thus, the pure complementary energy principle, first proposed in 1998 [3], leads to the fol¬ 
lowing canonical dual variational problem 

(V d ): sta|n d (C) = ^G d (C)dfi| CScS a |. (23) 

Since the canonical dual variable ( is a scalar-valued function, the criticality condition for 
this variational problem leads to a so-called canonical dual algebraic equation (see 0 ): 

4C 2 VU*(C) = T 2 (x) Vx £ B. (24) 


Note that VU*(C) :£*—>■ £ a is also one-to-one and onto, this equation has at least one 
solution for any given r 2 = tr(T 7 T) > 0 and Q = 0 only if r = 0. Therefore, although there 
is an inverse term in G d (Q, this canonical dual function is well-defined on S a . Due to the 
nonlinearity, the solution to (24) may not be unique nmm- By the pure complementary 
energy principle proposed by Gao in 1999 (see 0 ), we have 
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Theorem 2 (Complementary-Dual Principle) For any given T E T a , the following state¬ 
ments are equivalent: 

1) (xv) is a stationary point ofE(x,C)l 

2) x is a stationary solution to ( V); 

3) g is a stationary solution to (V d ). 

Moreover, we have 

n(x) = H(*,?) = n d (f) (25) 

Proof. For any given T E 7 " a , the stationary condition of 3(x, C) leads to the canonical 
equilibrium equations 


A(F(*)) = W*(?), 
2m*) = T E T a 


(26) 

(27) 


By the canonical duality, ( [26]) is equivalent to 7 = VP(£) with £ = A(V%). Thus, x must 
be a stationary solution to (V)t and also a stationary solution to CP) due to Theorem [I] 

By solving (27) we have F(x) = ^T. Substituting this into (26) leads to the canonical 
dual equation (24). Thus, 7 is a stationary solution to (V d ). 


The equivalence and the equation (25) can be proved by 


sta{n T (Vx)| X E Xc} = sta{3(x,C)| (x, 0 e T a x 4*} = sta{II d (C)l C G <S a } 

and Theorem [l] 


□ 


Theorem 3 (Pure Complementary Energy Principle) For any given nontrivial t / 0 
and x G X a such that T E T a / 0, (24) has at least one solution (k 7= 0, the deformation 


gradient defined by F& = VXk = C k 'T a critical point ofU(x) and II(x fc ) = n d (£fc)- 


Moreover, ifVx (C X T) = 0, then the deformation vector defined by 


Xk( x ) = \ l C fc X T • dx 

Z «/x 0 


(28) 


along any path from xo £ S x to x. £ B is a solution to (BVP )t in the sense that it satisfies 
both equilibrium equation and boundary conditions in m 


Proof. By the canonical duality relations in (18) we know that (k = X7V*((k) > 0. Thus, 
for a given nontrivial t(x), there exists a nontrivial t 2 (x) = tr(T 2 T) in B such that the 
canonical dual algebraic equation (24) have at least one nontrivial solution Cfc( x ) hr £>. 


Since the critical point Qk is a solution to (24), we have 


1 


f k = tr(Fj£Ffc) = s Crtr(T J T) = W*(C*0 


F A: = 2 C lT 


(29) 


subjected to any given rotation field R(x) E 77. By the fact that the canonical dual solution 


(k defined by (24) is independent of the rotation field, the canonical duality leads to 


G d (( z ) = 3(F*, C fc ) = F(A(Ffc)) - tr(F^T) = G(F k ). 

This shows II(x fc ) = B d ((k)- 

To prove Xk defined by (28) is a solution to (BVP) t, we simply substitute Vx fc = F^ = 

1/—lr 


X T into (BVP )t to have all necessary equilibrium conditions satisfied. Therefore, Xk 
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defined by (28) is a solution to ( BVP ): 


□ 


This pure complementary energy principle shows that by the canonical dual transfor¬ 
mation, the fully nonlinear partial differential equation in (BVP )t can be converted to an 


algebraic equation (24), which can be solved to obtain a complete set of solutions (see BIB!). 

Since S a is nonconvex, in order to identify global and local optimal solutions, we need the 
following convex subsets 

5+ = {C e 5 a | C > 0}, S~ = {C e S a \ ( < 0}. (30) 

Then by the canonical duality-triality theory developed in [5J we have the following theorem. 


Theorem 4 Suppose that V : < S a —>■ M is convex and for a given T E T a such that {£*,} is a 
solution set to ( 24 ), F& = |C A T 1 T, andXk defined by (28), we have the following statements. 


1. If( k £S+, then\7 2 W(F k ) >- 0 and x k a global minimal solution to (V) 

2. If Ck E Sf an d V 2 IT(F k ) >- 0, then x k a local minimal solution to (V). 

3. If C k E S~ and \7 2 W(F k ) ~< 0, then x k a local maximal solution to (V). 


If {Ck} C <S+, then {Xfc} is a convex set. The solution of (V) is unique if {Ck} C <S+. 


Proof. By using chain rule for W(F) = F(£(F)) we have VVF(F) = 2F[VF(£)] = 2£F, and 

V 2 1T(F) = 2(1 <8> I + 4/j(£)F ® F, (31) 

where I is an identity tensor in M 3x3 , h(C) = V 2 F(£) > 0 due to the convexity of V on £ a . 
Therefore, V 2 IT(F fc ) >- 0 if Ck E S+. 

To prove x k is a global minimizer of ( (P ), we follow Gao and Strang’s work in 1989 j 12j . 
By the convexity of V (f ) on its convex domain £ a , we have 

v(C) - v(c k ) >(C- C k )Ck V£, Ck € £ a , Ck = w(Ck)- (32) 


For any given variation Sx , we let x = Xk + $X- Then we have 

A(Vx) = tr[(Vx) T (V X )] = A(V Xk ) + A t (Vx fe )(V<Sx) - A C (V<5*), 


(33) 


where A t (F)5F = 2tr[F T (5F)] and A c (Sx) = -A(<5*). Clearly, A(F) = A t (F)F + A C (F). 
Then combining the inequality (32) and (33), we have 


n(x)-n(x fe )> / 2C k tr[(Vx k y(V8x)}dB - / 5x-tdS+ CM(V X f (^x)]dB 

JB JS t JB 

= [ [2Cfc(Vx fc ) - T] : (V6x)dB + G ap (5 X , Ck) V*, € X c (34) 

JB 

for any given T£7^, where 


G ap (x, 0= [ ~A c (V X )CdB= [ Ctr[(V*) T (Vx)]d B (35) 

JB JB 

is the complementary gap function introduced by Gao and Strang in m- If x k is a critical 
point of Il(x), then we have 


[ (2(V Xk )Ck - T] : (V8 X )dB 

JB 


o V<$x E T c , VT E T a 
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Thus, we have IT(x) - n(x fe ) > G ap (8x,Ck) > 0 V5x G X c if ( k G 5+. This shows that Xk is 
a global minimizer of (V). 


To prove the local extremality, we replace F& in (31) by = ^C k such that 


G(Cfc) = V 2 W(F k ) = 2( k I <8 I + C 2 h{i k ) T ® T, 


(36) 


where f k = W*(Ck)- Clearly, for a given T G T a such that Ck G S~, the Hessian V 2 IT(Ffc) 
could be either positive or negative definite. The total potential n(x fc ) is locally convex if 
the Legendre condition V 2 W (\7 Xk) G 0 holds, locally concave if V 2 IT(V;x fc ) -< 0. Since Xk 
is a global minimizer when Ck G Sfi, therefore, for Ck G S ~, the stationary solution x k is a 
local minimizer if V 2 lT(Vx fc ) >- 0 and, by the triality theory^, HT] , x k is the biggest local 
maximizer if V 2 IT(Vx fc ) -< 0. 

If {Cfc} C Sfi, then all the solutions {x fc } are global minimizers and form a convex set. 
Since n d (£) is strictly concave on the open convex set «S+, the condition {£*.} C <S+ implies 
the unique solution of (|24|). In this case, Problems (V)t has at most one solution. □ 


Theorem 5 (Triality Theory) For any given T G 7^ / 0, let Ck be a critical point of (V° 
the vector Xk be defined by (28), and X c 
If ( k G <S+, then 


x S 0 C X c x S a a neighborhood of (x k , C k)- 


n (Xfc) = min n(x) 

If Ck G Sf and G(Ck) >~ 0 , then 

n (Xfc) = min n(x) 

A 

If Ck G and G(Ck) ~< 0, then 


n (Xfc) 


max n(y) 
XGT, 


maxn d (C) = n d (C fc )- 

(37) 

CS-S+ 


mmn d (C) = n d (Cfc). 

(38) 

C eSo 


maxn d (C) = U d (C k ). 

C 

(39) 


This theorem shows that for convex canonical function V, the triality theory can be used 
to identify both global and local extremum solutions to the variational problem (V) and 
the nonconvex minimum variational problem (’P)t is canonically equivalent to the following 
concave maximization problem over an open convex set Sfi. i.e. 

(iPtt) T : max jn rf (C) = J G d (C)dB | C G 5+1 , (40) 

which is much easier to solve than directly for obtaining global optimal solution of (V). 


3 Generalized Quasiconvexity, G-Ellipticity, and Uniqueness 

Ellipticity is a classical concept originally from linear partial differential systems, where the 
deformation is a scalar-valued function x '■ £> —> M and stored energy is a quadratic function 
W(~f) = 7 of 7 = Vy G M 3 . The linear operator 

L(x] = -V • [H(Vx)] = -[hijXdli 

2 The neighborhood X 0 of Xk * n the canonical duality theory means that Xk i s the only one critical point 

of n(x) on X 0 (see [5]). 
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is called elliptic if H = {hij} is positive definite. In this case, the function G( 7 ) = W ( 7 )— 7 T t 
is convex and its level set is an ellipse for any given r £ M 3 . This concept has been extended 
to nonlinear analysis. The fully nonlinear partial differential equation in ( BVP ) ([h]) is called 
elliptic if the following Legendre-Hadamard (LH) condition holds 

(a (g) a) : V 2 IV(F) : (77 <g) 77 ) > 0 Va, 77 £ M 3 , VF £ T a . (41) 


The (BVP) is called strong elliptic if the inequality holds strictly. In this case, (BVP) has 
at most one solution. In vector space, the LH condition is equivalent to Legendre condition 
V 2 IL(7) 7 0 V7 £ M n . 

Clearly, the LH condition is only a sufficient condition for local minimizer of the variational 
problem ( V). In order to identify ellipticity, one must to check LH condition for all local 
solutions, which is impossible for general fully nonlinear problems. Also, the traditional 
ellipticity definition depends only on the stored energy IT(F) regardless of the linear term in 
G(F) = VF(F) — tr (F 7 T). This definition works only for convex systems since the linear term 
tr (F 7 T) can’t change the convexity of G(F). But this is not true for nonconvex systems. To 
see this, let us consider the St. Venant-Kirchhoff material 


W(F) = : H : E, E = ^[(F) T I 


F) - II 


(42) 


where I is a unit tensor in M 3x3 . Clearly, this function is not even rank-one convex. A special 
case of this model in M n is the well-known double-well potential W( 7 ) = ^(^| 7| 2 — l) 2 . 
If we let £ = A( 7 ) = AI 71 2 — 1 be an objective measure, we have the canonical function 
V(£) = |£ 2 . In this case, the canonical dual algebraic equation (24) is a cubic equation 
(see m 2 <^ 2 (C + 1) = t 2 1 which has at most three real solutions {£/(x)} at each x £ B 
satisfying £1 > 0 > C 2 > C 3 - It was proved in [5] (Theorem 3.4.4, page 133) that for a given 
force t(x), if r 2 (x) > 8/27 Vx £ B C M, then (BVP) t has only one solution on B. If 
t 2 (x) < 8/27 Vx £ B s C B, then (BVP) t has three solutions {xfc( x )} at each x £ B s , i.e. 
n(x) is nonconvex on B s ■ It was shown by Gao and Ogden that these solutions are nonsmooth 
if r(x) changes its sign on B s [El- 

Analytical solutions for general 3-D finite deformation problem (V) were first proposed 
by Gao in 1998-1999 [3j 2f. It is proved recently [8] that for St Venant-Kirchhoff material, 
the problem (V) could have 24 critical solutions at each material point x £ B, but only one 
global minimizer. The solution is unique if the external force is sufficiently large. 

For a given function G : T a —> M, its level set and sub-level set of height a £ M are 
defined, respectively, as the following 


C a (G) = {F £ Ta | G(F) = a}, C a {G) = {F £ T a \ G( F) < a}, a £ 


(43) 


The geometrical explanation for ellipticity and Theorem [4] is illustrated by Fig. [lj which 
shows that the nonconvex function G( 7) = ^(^ I7I 2 — l) 2 — 7 t t depends sensitively on the 
external force r £ M 2 . If |r| is big enough, G( 7 ) has only one minimizer and its level set is 
an ellipse (Fig. [T] (b)). Otherwise, G( 7 ) has multiple local minimizers and its level set is not 
an ellipse. For r = 0, it is well-known Mexican-hat in theoretical physics (Fig. HIM). 


Fig. a shows that although G( 7 ) has only one global minimizer for certain given r, the 
function is still nonconvex. Such a function is called quasiconvex in the context of global 
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Figure 1: Graphs and level sets of G(x) for r = 0 (left) and r / 0 (right) 


optimization. In order to distinguish this type of functions with Morry’s quasiconvexity in 
nonconvex analysis, a generalized definition in a tensor space J~ a C W nxn could be convenient. 

Definition 1 (G-Quasiconvexity) A function G : T a C M mxn — y M is called G-quasiconvex 
if its domain T a is convex and 

G{6 F + (1 - 0)T) < max{G(F), G(T)} VF, T € T a , V6 E [0,1], (44) 

It is called strictly G-quasiconvex if the inequality holds strictly. 

Moreover, we may need a definition of generalize ellipticity for nonconvex systems. 

Definition 2 (G-EIlipticity) For a given function G : T a —> M and a E M, its level set 
C a {G) is said to be a G-ellipse if it is a closed, simply connected set. For a given t such that 
T E T a , the {BVP) is said to be G-elliptic if the total potential function G(F) is G-quasiconvex 
on T a - {BVP) is strongly G-elliptic if G{ F) is strictly G-quasiconvex. 

Lemma 1 For a given function G : T a C M mxn — y M ; 

G(F) is G-quasiconvex £„(G) is convex C a {G) is a G-ellipse Vq £l . (45) 

G(F) is convex=$■ is rank-one convex=$■ G-quasiconvex {BVP) is G-elliptic. (46) 

This statement shows an important fact in nonconvex systems, i.e. the total number of 
solutions to a nonlinear equation depends not only on the stored energy, but also (mainly) on 
the external force field. The nonlinear partial differential equation in {BVP) is elliptic only 
if it is G-elliptic. {BVP) has at most one solution if G(F) is strictly G-quasiconvex on T a . 

Remark 3 (Existence and Uniqueness) Suppose that the canonical function V : £ a —> M 
is convex, then VU*(£) > 0 is a monotonic operator on £*. If for any given t : St —> M 3 such 
that T G T a / 0 and r 2 (x) = tr( T 2 T) / 0 Vx G B. then the nonconvex variational problem 
{V) has at least one nontrivial solution a.e. in B. It has a unique nontrivial solution if there 
exists a constant t c such that t 2 (x) = fr(T 3 T) > t 2 Vx G B. 
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In global optimization, the most simple quadratic integer programming problem 

0 P)i : min jll(x) = Jx T Qx - x T t | x = {xi} n G {0, l} n C R n 

could have up to 2 n local minimizers, which can’t be solved directly by traditional determin¬ 
istic methods in polynomial time due to the indefinite matrix Q and the integer constraint. 
Such a nonconvex discrete optimization problem is considered as NP-hard in computer sci¬ 
ence. However, by using canonical transformation £ = A(x) = {xi(xi — 1)} G M n , the 
canonical dual of this discrete problem is a concave maximization over a convex set in contin¬ 
uous space m- It was proved in | 6 j that there exists a positive vector t = {rj} n >06 M n , 
if {|fi| < n} n , then 5+ / 0 and (V)i is not NP-hard. The decision variable is simply 
{xi} = {0 if ti < —Ti, 1 if ti > Ti} (Theorem 8 , [ 6 ]). Thus, the canonical duality theory can 
be used to identify NP-hard problems mi- 


4 Applications in Anti-Plane Shear Deformation 

Now let us consider a special case that the homogeneous elastic body B C M 3 is a cylinder 
with generators parallel to the e 3 axis and with cross section a sufficiently nice region Q C M 2 
in the ei x e 2 plane. The so-called anti-plane shear deformation is defined by (see PI) 

x(x) = {x\, x 2 , x 3 + u( Xl ,x 2 )} :!1cM 2 4M 3 , (47) 

where (xi,x 2 ,x 3 ) are cylindrical coordinates in the reference configuration B relative to a 
cylindrical basis {e,}, z = 1,2,3. On T x C <912, the homogenous boundary condition is given 
u(x a ) = 0 \/x Q G r x , a = 1,2. On the remaining boundary T t = d^l 0 T x , the cylinder is 
subjected to the shear force 

t(x) = t(x)e 3 Vx G T f , 

where t : Tt —> M is a prescribed function. For this anti-plane shear deformation we have 

/ 1 0 0\ /1 + Uj ^,i\ 

F = Vx=|0 10, C = F 1 F=| u t iu }2 1 + u 2 2 u y2 I , (48) 

\U, 1 U } 2 1 / \ U 1 U 2 1 / 

where u tCt represents du/dx a for a = 1,2. By the notation |Vu| 2 = u 2 ^ + u 2 2 , we have 

Ii(C) = / 2 (C) = 3 + |Vu| 2 , 7 3 (C) = 1, (49) 


Clearly, both F and /i(C) depend only on the shear strain 7 = Vu = {u iQ }, therefore, the 
strain energy can be equivalently written in the forms of 


W( F( 7 )) = V(Ch)) = Wh) 


(50) 


where W ( 7 ) is a real-valued function. 

The fact detF = 1 shows that the anti-plane shear state (47) is an isochoric deformation. 
Therefore, the kinetically admissible displacement space X c can be simply replaced by a 
convex set 

U c = (u(x) G yV 1 ’ 1 (P;M)| rt(x) = 0 Vx = {x a } G T x }. (51) 
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Thus, in terms of £ = A(j) = I\ — 3 = |^y| 2 and W(F( 7 )) = for any given 


t £ Ta = {r E C 1 [f2;R 2 ]| V-r = 0 in 0, n -r = t on T t } 


Problem (V)t f° r the anti-plane shear deformation (47) has the following form 


(V) s : min } II(tt) = / G(Vu)df2 | u E U c } , G{~f) = F(A( 7 )) — 'y 1 r (52) 

In 


Under certain regularity conditions, the associated mixed boundary value problem is 


(BVP) S : 


( V • (2CVit) = 0 in O, 

} n • (2£Vn) = t on r ( , u = 0 on T x 


(53) 


where n = {n a } £ M 2 is a unit vector norm to <912, and £ = VV(£), £ = |Vu| 2 . 

If T x = <912, then ( BVP) S is a Dirichlet boundary value problem, which has only trivial 
solution due to zero input. For Neumann boundary value problem T^ = <9f2, the external 
force field must be such that 


/ t(x)dr = 0 

JT t 


for overall force equilibrium. In this case, if x is a solution to ( BVP ) S , then x = X + c is 
also a solution for any vector c E M 3 since the cylinder is not fixed. Therefore, the mixed 
boundary value problem ( BVP) S is necessary for anti-plane shear deformation to have a 
unique solution. 

By the fact that the only unknown u is a scalar-valued function, anti-plane shear defor¬ 
mations are one of the simplest classes of deformations that solids can undergo m- Indeed, 
if U(£) is a canonical function on £ a = {£ E L P (I2)| £(x) > 0 Vx E 12} and for any given 
r E Ta such that r = |t|, the canonical dual problem has a very simple form 


C r d )s : sta|n d (C) = ^ 


-no-^rv 


df2 | £ E «S D 


(54) 


Since A(u) = |Vu| , the canonical dual algebraic equation (24) for this problem is 


4C 2 VU*(C) = r 2 (x), VxEfi. 


(55) 


Corollary 1 For any given non-trivial shear force t(x) ^4 0 on Ft such that r E 7^ / 0, the 
canonical dual problem (' P d ) s has at least one non-trivial solution Ck- I/V x Ck T = 0, the 
scale-valued function 

1 /’ x 

«fc(x) = 9 / CV • dx (56) 

along any path from xo E F v fo x E £2 is a critical point o/II(u) and Il(ufc) = II rf (£fc). 

If Ck £ ; then Uk is a global minimizer of (' V) s . 

If Ck S S~ and G(£fc) 7 0, then Uk is a local minimizer of (' P) s . 

If Ck S S~ and G (Ck) P 0, then Uk is a local maximizer of (V) s . 
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Example. Applications of the canonical duality theory to general anti-plane shear problems 
have been demonstrated for solving convex exponential and nonconvex polynomial stored en¬ 
ergies recently in j7J. In this paper, the following generalized neo-Hookean model is considered 

V(0 = ci (h - 3) + c 2 {h - 3) log(/i - 3) (57) 

where ci, c 2 are positive material constants. Clearly, E(£) is convex in £ = I\ — 3, but 

W{l) = V(h{ 7 )) = ci|7| 2 + c 2 |7| 2 log|7| 2 
is a double-well function of the shear strain 7 = Vrt (see Fig. [ 2 ]). 




(a) 


(b) 


Figure 2: Graphs of W( 7 ) (a) and its derivative (b) (ci = c 2 = 1) 

It is easy to check 

C = W(0 = Cl + c 2 (log£ + 1) : S a ->• El = L q (Q) 

is one-to-one and onto, where q is a dual number of p > 1, i.e. 1/p + 1/q = 1. The 
complementary energy can be obtained easily 

V*(C) = sta{£C - ^(£)| C e £a} = c 2 exp[c^ 1 (C - ci) - 1] 

In this case, the canonical dual algebraic equation is 

"C - c i 


C 2 exp 


- 1 


c 2 


= r 2 (x) Vx G 12. 


(58) 


Let h 2 (Q = C/ 2 exp[(^ —ci)/c 2 —1] be the left hand side function in the canonical dual algebraic 
equation (58). By solving h'((, c ) = 0 we known that at Ce = —2c 2 , h(C) has a local maximum 

hmax(Ce) = V = 2c 2 y / exp[-3 - ci/c 2 ]. 

From the graphs of the canonical dual algebraic curve h(() given in Fig. [3]we can see that 
the canonical dual algebraic equation ( |58[ ) may have at most three real solutions in the order 
Of Cl > 0 > C2 > Cs depending on r = |r(x)|, x £ (see Fig. §r). The equation ( [58]) has a 
unique solution if r > rj. In this case, the total strain grand G(y) is strictly G-quasiconvex 
(see Fig. [i|. Fig [ 5 ] shows the graphs of G(y) and its canonical dual G d ( C) for r < p. In this 
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case, the function G{ 7 ) is nonconvex and has three critical points. The triality theory holds 
for G(j) and its canoncial dual G d (() 

G( 71 ) = min G( 7 ) = max G d (() = G d (( 1 ). 


G( 72) = minG(7) = min G d (() = G d (( 2 ). 

l&Qo C>-2C2 

G( 73) = maxG(7) = max G d (C) = G d (C 3 ), 
7600 C<-2c 2 

where is a neighborhood of 7$ (z = 1,2). 



Figure 3: Dual algebraic curve h(Q (ci = C 2 = 1) 




Figure 4: Graphs of G-quasiconvex G( 7 ) (ci = C 2 = 1) 


5 Conclusions 

In summary, the following conclusions can be obtained. 

1. The pure complementary energy principle and canonical duality-triality the¬ 
ory developed in j5] are useful for solving general nonlinear boundary value 
problems in nonlinear elasticity. 
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(b) Graph of G d (Q 


Figure 5: Graphs of G{ 7) and G d {Q for t < r/ (c\ = C2 = 1 ) 


2. Both convexity of the total potential and ellipticity condition of the associated 
fully nonlinear boundary value problem depend not only on the stored energy 
function, but also sensitively on the external force field. 

3. The Legendre-Hadamard condition is only a necessary ellipticity condition for 
convex systems. The triality theory provides a sufficient condition to identify 
both global and local extremum solutions for nonconvex problems. 

These conclusions are naturally included in the canonical duality-triality theory developed 
by the author and his co-workers during the last 25 years [5]. Extensive applications have 
been given in multidisciplinary fields of biology, chaotic dynamics, computational mechan¬ 
ics, information theory, phase transitions, post-buckling, operations research, industrial and 
systems engineering, etc. (see recent review article ED- 
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